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We investigate numerically the dynamics of a single polymer in a linear shear flow. The efi'ects of 
thermal fluctuations and randomly fluctuating velocity gradients are both analyzed. Angular, elon- 
gation and tumbUng time statistics are measured numerically. We perform analytical calculations 
and numerical simulations for a linear single-dumbbell polymer model comparing the results with 
previous theoretical and experimental studies. For thermally driven polymers the balance between 
' relaxation and thermal fluctuations plays a fundamental role, whereas for random velocity gradients 

' the ratio between the intensity of the random part and the mean shear is the most relevant quantity. 

I In the low-noise limit, many universal aspects of the motion of a polymer in a shear flow can be 

04 . understood in this simplified framework. 
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PACS numbers: 83.80.Rs,47.27.Nz,36.20.- 



I. INTRODUCTION 

, Thanks to recent irnprovements in experimental techniques it is nowadays possible to follow the motion of individual 

■ molecules in solvents 0,|3;I1;SISIESIE0I- The characterization of polymer dynamics at the level of a sin gle m olecule 
• I is a first step towards the understanding of mechanical interactions between biomolecules (see e.g. 'To',^^ 0, 0, 

■ 0j ^8 JJ] ) ; of the fundamental rheology of polymer solutions, and of the viscoelastic properties of more complex 
flows (see for example and references therein for elastic turbulence). 

Recently measurements of elongation and orientation with respect to simple external flows have been performed 0, |^ 
in order to analyze how the conformation of a single molecule can be modified by an external field. 
" ' , As the number of degrees of freedom needed to fully describe a macromolecule is extremely high it is necessary to 
' formulate theoretical models, to verify them with numerical simulations, to simplify the problem , and to understand 
which observables play a key role. 

Large numbers of papers have been written on the subject of single polymer dynamics in shear flows from experimental, 
^\ _ numerical and theoretical viewpoint by Chu, Larson, Shaqfeh and their respective collaborators. 

i A polymer molecule in a plane, linear, steady shear flow is oriented in the flow direction by the velocity field. As the 
' rotational motion of the polymer is determined only by the velocity difference in the space points, when it is aligned 
. along the shear direction the external flow effect becomes negligible, and the thermal noise is the most relevant 
P5 ■ external force. Thermal fluctuations can skew the polymer to regions where the external flow is relevant again. In 
these cases the shear flow can induce a fast rotation and align the polymer again along the (reverse) flow direction, 
^ i.e. the polymer tumbles H S 0, II II EH ■ 

^ [ This phenomenon can happen via several conformational pathways due to the complexity of the motion of a polymer 
•T-j . molecule (see for example ]1,I3), and can be fully described only by taking into account all the degrees of freedom of 
^\ ' the polymer. Unfortunately in the framework of these complex polymer models it is very difficult to obtain analytical 

^ I results. 

- - • The main goal of our paper is to test numerically the predictions of recent theoretical and experimental studies 0, HO] 
in a framework in which analytical results can be obtained 1^. and to analyze the statistics of the tumbling times, 
i.e. the time between two subsequent flips of the polymer [3ll7ll2l|. The simplest model that reproduces qualitatively 
the behavior of polymers is the dumbbell model |23| . This model allows one to carry out some analytical calculations 
in the case which we analyze ^3 113, 113, 13 > is veryeasy and fast to simulate numerically (see sec. Inland ^25j) and 
reproduces qualitatively recent experimental results 0, ■ 

The paper is organized as follows: in sec. ^ the evolution equation of the polymer and the numerical methods are 
briefly explained. Sec. lIIIl is devoted to the analysis of the statistics of thermal fluctuations of a flexible polymer placed 
into a linear shear flow. In our work we present the analysis of the stationary distribution of the polymer end-to-end 
vector and we study the distribution function of the polymer tumbling time, which can be measured experimentally. 
In sec. liVl we study the angular dynamics of strongly elongated polymers, for which the size fluctuations are negligible. 
Finally in sec. we study the elongation statistical properties of the end-to-end vector in a random flow with a large 
mean shear. 
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II. BASIC RELATIONS AND NUMERICAL ANALYSIS 

We wish to analyze the behavior of a polymer in a generic simple shear flow experiencing the Langevin force |l9j |. 
In general there are two effects of the velocity field on a polymer molecule: the Lagrangian advection of the polymer 
and the elongation/relaxation dynamics due to velocity gradients. In all the cases discussed in this paper we disregard 
the Lagrangian dynamics by staying in the reference frame of the polymer center of mass. For the internal degrees 
of freedom of the polymer we use the simple dumbbell model ^22^, leaving the analysis of more realistic models for 
future studies |2^. In this case the basic equation describing the evolution of the polymer end-to-end vector R has 
the following form: 



= a,,Rj -jR^ + \l (1) 

where ct^ = djVi is the velocity gradient matrix, 7 is the polymer relaxation rate, ^ is the thermal noise term, which 
has white-noise statistics: (^i(i)Cj(i')) — 5ij5{t — t'), and Rq is the equilibrium length in the absence of an external 
field. 

For incompressible flows there is no ambiguity in the discretization of the stochastic differential equation and 
no additional contact term must be taken into account, so we resort to the Euler-Ito scheme without any loss of 
generality. 

When the gradient a contains only a steady linear shear contribution we can write down the formal solution of 
as: 

R,{t) = i?,(0)e-''* + e--^* /* e''*'s,jRj{t')5,,6jydt'+ 



-7e-^* e-"' R.it')dt' + ^^e"-^* e"'Ut')dt'. (2) 

This equation can be discretized, and the terms containing the thermal noise can be rewritten as new gaussian 
variables with amplitudes that can be computed directly, so that the flnal solution reads: 



J_(l _ e-27At) _ iAie-27At _ At2e-27At 
272 7 



(3) 



where {ilirfj) = 5ij6^^ . The subscripts stand for the cartesian coordinates while the superscripts stand for the 
discretized time. As expected in the limit Ai — )■ there is no contribution of rjy in the first equation and the 
amplitudes are the same as in eq. (^. 

In the case of a random velocity gradient, Wij , we have to generate the variables w^j such that they have the prescribed 
correlation function: 

{(T..j{t)aki{t')) = D5{t - t'){AS,k5ji - Sa&kj - > (5) 

where again the Dirac delta function is substituted by a Kronecker symbol. In this case the modulus of the vector 
R grows indefinitely so that we can normalize it at each time step. In order to describe the elongation properties we 
can compute the maximum Lyapunov exponent A and the corresponding finite time Lyapunov exponent At |27| . 

III. THERMALLY DRIVEN POLYMERS 

In this section we will examine the case of a linear steady shear flow in the plane XY: Uij = sSixSjy, where s is the 
shear rate. Eq. |^ has an explicit solution of the following form: 

R,{t) = exp{-jt)W,j{t)Rj{0) + f dt'e^v [-7(t - t')] W,,{t - t%{t'), (6) 
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FIG. 1: Schematic figure explaining polymer orientation geometry. 



where W{t) = exp(icr). At large times the initial polymer elongation is forgotten and after averaging over the thermal 
fluctuations ^ one can easily obtain the following distribution function: 

P{R) = (2^)-3/2(det/)-i/2exp (7) 
^ = ^ /r eM-'2lt)W{t)W^it) . (8) 

The probability density function (PDF) Q is valid for any velocity gradient. In the particular case of a steady shear 
flow in the XY plane, as shown in Fig. ^ the matrices can be found explicitly: 

/ I St 0\ n2 /I + SV272 s/27 0\ 

Wit) =010, / = ^ s/27 10, (9) 

yooiy •^y o 01/ 

where the axis are sorted in the order X,Y,Z. One can see that for large Weissenberg numbers Wi = s/7 ^ 1, the 
mean polymer elongation in the X direction is much larger than in transversal directions. 



A. Elongation PDF 



The PDF of elongation can be easily computed from (|7|): the elongation in the y and z direction are independent 

of Wi and the two marginal PDFs are gaussian with zero mean value and variance (Rl) = {Rfi = ^, while in the 

mean flow direction the variance is (i?^) = -^(1 + ^Wi^). 

The distribution function of R cannot be obtained analytically for an arbitrary value of Wi, but one can study its 
behavior in two limiting cases. For Wi <C 1 the effect of the shear flow can be neglected, and we have the simple 
thermally driven polymer with linear relaxation force. The distribution function of its elongation has a simple gaussian 
form I2I: 



2i?2 



2Rl 



(10) 



In the opposite case Wi 3> 1 the system is strongly anisotropic and the main contribution to the polymer elongation 
comes from the X component. This fact allows one to obtain the right tail of the elongation PDF. For R^ Rq one 
has 



P{R) = 



12 



1 



7r(2 + Wi^) -Ro 



exp 



(2 + Wi^)i?2 



(11) 



B. Angular PDF 



Next we are interested in the orientation of the polymer. The angular distribution function will allow us to show 
that for large Wi the polymer spends most of the time aligned to the X axis, while for Wi <C 1 the orientation 
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FIG. 2: The PDF of the modulus of the end-to-end vector R for different Wi numbers. While for Wi = the PDF is 
concentrated around the radius of equilibrium, for higher values of Wi it is broader and the high elongation configurations 
become more and more probable. 
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(a) The analytic PDF is the numerical integration of 
I12i . while the dots come from the numerical 
simulations. 



(b) The amplitude of the PDF's peak i 
of Wi 



as a function 



FIG. 3; The PDF of the angle 



distribution is almost isotropic. In order to parameterize the orientation of the polymer we introduce the angles (j) and 
as shown in Fig. ^ The angle (p represents the deviation of the polymer end-to-end vector from the X axis in the 
shear velocity plane XY while the angle 9 gives us the polymer deviation in the transversal direction Z . Switching to 
spherical coordinates we have Rx = Rcos9 cos<j), Ry — i? cos 6* sin 0, Rz ~ Rsm9. After integrating over the polymer 
elongation in eq.Q one immediately obtains the angular PDF: 



cos 9 



1 



cos^ 6 
4+W? 



[Wi^ cos(20) + 2Wi sin(20)] | 



3/2 



(12) 



The calcul ation of the two marginal PDFs is not possible in general, and we have to integrate eq. (|12|l numerically 

(Fig.|gi3. 

As one can see in Fig.01 the PDF of 9 decays as a power law (asymptotically 1^(9) ^ 9 ^ in the region (pt -€1 9 <^ 1), 
and is symmetric about 9 = 0. The higher the value of Wi, the wider is the power law region. 

In principle, if the polymer tumbles in the shear plane, we should consider a time-dependent PDF of on a unbounded 
domain. However if we consider (f) to be between and 27r one arrives at a stationary PDF, peaked in the regions 
(j) = kn 

V{4>) ~ E 



(j)t,k = 0,1,2 (Fig. 3(a) I, where 

,-2 / 



as shown in Fig. 3(b) In the region where (j) ^ cpt we have 
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FIG. 4: The PDF of the angle 9 is symmetric with respect to 9 = Q and can be described by a power law relationship over a 
wide range of values of 9. 




FIG. 5: The two different tumbling time PDFs. 



The angular PDFs tell us that any fluctuations in 6 are not relevant to the tumbling dynamics, and that for the angle 
there is no relevant scale. The fact that increasing Wi the power law region becomes wider and wider means that 
when the shear is much stronger than the relaxation the polymer is nearly aligned with the X axis. The PDF of the 
angle expresses the fact that the polymer spends most of the time in the vicinity of = 0, = tt, (/> = 27r. The 
polymer is aligned in the X direction as a result of the shear flow, it fluctuates for a certain time in the region </) </)(, 
where the shape effects are strong. When thermally activated it goes beyond the region </> ^ </>( and the shear becomes 
more important and induces a fast rotation of A(/> ~ tt or multiples thereof, that is the polymer tumbles. 



C. Tumbling time distributions 



In the previous two sections we have studied the stationary distribution functions which can be measured exper- 
imentally by averaging the polymer elongation and orientation over large time periods. However, as described in 
the introduction and in |20l l23j | the dynamics of the polymer is non stationary, due to continuous tumbling. The 
natural question which arises is whether there are some quantities, which would allow experimental observation and 
quantitative description of the tumbling process. In the recent papers by Chu et al. the different time correlation 
functions of the polymer elongation were studied, and it has been shown that they can have different forms in the 
presence of shear. Here we introduce the tumbling time, which can be measured experimentally and used to 
describe the tumbling process. As described above the most striking difference between the tumbling dynamics of 
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the polymer in a shear flow and spinning in a rotational flow is the a-periodicity of the tumbling process. Due to 
the stochastic nature of tumbling, the tumbling time, i.e. time between subsequent flips of the polymer, is a random 
variable with relatively large fluctuations. Our task is now to study the distribution function of the tumbling time 
and to give some theoretical and numerical predictions of its dependence on the Weissenberg number. 

In the case of the dumbbell model, we can define tumbling as a flip of the polymer induced by the thermal noise. 
When the polymer is in an unstable equilibrium configuration cj) — kir^k — 0,1,2 the thermal noise can bring the 
polymer out of the region (f) ^ (pt, and this induces a fast rotation which takes a time of the order of s^^. We define 
the tumbling time as the interval between two subsequent fiips. This is typically the time between subsequent 



crossings of the lines = (fc + 1/2)tt. For Wi 3> 1, as shown in the previous section (see Fig. 3(a) I, the polymer 
spends only a small fraction of the total time far from the shear direction, and therefore the exact value of the angular 
threshold defining the tumbling time is not very important. Indeed this time is made up of two contributions: the 
time spent in the region (j) ^ (pt and the duration of the rotation of order s^^. The latter can be neglected for large 
Wi. 

Experimental techniques do not always allow the polymer orientation angle (j) to be resolved and therefore it is 
convenient to introduce another definition of the tumbling time {rji). This quantity measures the interval between 
subsequent changes between stretched and coiled polymer conformations. In other words, we start measuring when 
the length of the polymer exceeds a certain threshold value and we stop when it again becomes smaller than this 
threshold value. 

While T0 is the most natural definition it can be difficult to measure experimentally. On the other hand, the PDF of 
Tn depends in a non universal way on the threshold, but experimental techniques allows it to be measured when the 
polymer is sufficiently stretched. 

In order to analyze the behavior of the tumbling time PDFs at large Weissenberg numbers we should analyze 
the dynamics of the polymer elongation projected onto the X-axis (Rx)- Indeed, the angle-based tumbling time 
measures the time intervals between subsequent events when = 0, while the elongation based tumbling time th 
corresponds to the time intervals between crossings of a threshold value Rx — Rth- In order to find how the shape 
of the tumbling time PDF changes with increasing Wi number, we can rcscalc the time and the polymer size in the 
following way: R^ — xRoWi, Ry — yRo, t = 7~^t, which leads to the following equations: 

X = —X + y (13) 
2 



^Cy (14) 
(C,(0)C,(r))=<5(r). (15) 

In eq. (O we have omitted the term corresponding to S,x because it is negligible in the large Weissenberg number 
regime compared to the other terms. Note that these equations do not depend on the Weissenberg number, and 
it enters the problem only through the threshold Rth- In our simulations we choose Rth = \/3(i?^) cx i?oWi, so 
that in the dimcnsionless variables we have xth = 1/V2 which does not depend on Weissenberg number either. 
From this analysis we conclude that the tumbling time PDF at large Weissenberg numbers approaches some universal 
form: its peak is positioned at characteristic time scale of order 7~^, and as 7 increases the peak moves towards 
the origin. The problem of finding the exact form of this function is equivalent to the problem of finding the 
PDF of persistence times of a non-Mar kovian random process x{t). This problem has recently attracted a lot of 
attention (see e.g. ref. H^), however no explicit solution is known for a general random process. Still it is possible 
to make some predictions on the PDF: the right tail of the PDF tj^, ^ 7^^ is exponential and has the form 
PijR) oc exp [—CfCiTRl, P{T(t,) (x exp [—c^jt^] because large tumbling times correspond to a large number 7r ^ 1 of 
unsuccessful attempts to cross the threshold. The correlation time of our stochastic process Rx{t) is of order 7"^, 
and therefore these attempts to cross the threshold are almost independent, and one should take the product of their 
probabilities, which leads to the above exponential laws. Both the PDFs of and are not well defined for ^ 0, 
Tji since they are sensitive to brownian noise discretization: it is possible to tumble very rapidly when the polymer 
is coiled, and this is why the probability of measuring a tumbling time much smaller than is not zero. Both the 
tumbling time definitions do not work very well in the case of small Wi experimentally because of the high resolution 
needed, and numerically because of the discretization procedure. 



IV. STRONGLY ELONGATED AND RIGID MOLECULES 



Another physical situation we are interested in is the dynamics of strongly elongated polymers in random fiows. 
In this model, described in detail in [20I I2M |2^ . the polymer is placed into a random fiow above the coil-stretch 
transition, where the effect of the thermal fluctuations can be neglected. In this case the orientational dynamics of 
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the polymer are decoupled from the evolution of the elongation, so that we can introduce the unit vector rii = Ri/R, 
obeying the following evolution equation: 

hi = nj{6ik - nink)VjVk- (16) 

The velocity gradient consists of a regular shear part (as in the previous section), and an isotropic incompressible 
random part. The polymer size in the experiments is always much smaller than the characteristic viscous scale of 
the velocity field which allows us to assume smoothness of the velocity field. We will also assume short-correlated 
velocity field which corresponds to the so called Kraichnan-Batchelor model, which has been extensively studied in 
recent years j29j . In the framework of this model the velocity gradient matrix aij — V^Wi in the Lagrangian frame is 
described by a gaussian process with the following pair-correlation function 30] : 

{(T^AtWl{t')) - D5{t - t'){A5,k5ji - 5^l5ko - . (17) 

where Uij is the gradient matrix of the random velocity component. Using the same notations as in the previous 
sections one arrives at the following dynamical equation for the angles: 

(/) = -ssin^ + ^0 , (18) 

= -s^^^^^sin0cos0 + Ce, (19) 

where and are zero mean random variables related to the fluctuating components of the velocity gradient. The 
statistics of both and can be obtained from the correlation function 117|l : 

{ieme{t'))=AD5{t-t') (20) 

mW)) = ^Q5{t-t'). (21) 

Note that the measure of configurations with 9 ~ 7r/2 is small, thus making the formal singularity in eq. I|21() 
not essential. From eq. H18|l it turns out that the polymer experiences constant a-periodic tumbling in the XY 
plane 12^12^1^. 



A. Stationary angular PDF 
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FIG. 6: The function U{x), where x = tanO/ sin(j}. 



From a dimensional analysis it can be shown that for s D the characteristic values of the angles 0, 9 will be of 
order {D/s^/^ <C 1. In this region one can set cos 9 = 1 in the expression l|21|) so that eq. (|18|l becomes completely 
independent of 9. Hence in this limit one derives: 

LO 

Psti'f) ^ Jj J d^exp 



AS 



8D^ ^' 24D 



(22) 
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(a) The agreement between the PDF and sin ^ <^ is 
good in the region </> 2> (/<t + kir, fc = 0, 1, 2 



(b) The behavior of (f>t as a function of D/s, with the 
dotted line shown as a guide only. 



FIG. 7: The PDF of the angle 



where uj is the mean rotation frequency of the polymer, which is determined from tlic normalization condition 
Jq Pst{4>)'^4> = 1 and is given by 



2U/3 



4-3i/6r(7/6)0F 



(23) 



The explicit form of the joint angular PDF is hard to compute analytically. However one can obtain the expression 
for the tails of the PDF 0, 6i > {D/s^/^ Elll: 



P( 



?7 (tan 6*/ sin 0) 
sin'^ 6 cos^ 6 



(24) 



where U{x) is an unknown function with an universal argument. Numerical simulations confirms this prediction, as 
shown in Fig. 

The two marginal PDFs obtained by numerical simulations show a behavior similar to that of thermally driven 
polymers. Also the phenomenology of the systems are similar, but now all the quantities depend on the ratio 
between the amplitude of the random velocity gradient and the mean shear. The PDF of the angle (j> has peaks at 
- (/)t + fcTT, A: = 0, 1,2, where c^t ~ {^Y^^ and the PDF of 6 has an algebraic core ViO) - 6^^ for (-|) < 1 (see 
Fig. 13). 

In the stochastic region |(/)| ^ (pt the tails of the 0-angle PDF P{d) are also algebraic P(0 = 0,9) oc 9^^, but the 
exponent y is non- universal and depends on the statistical properties of the random velocity gradient (see Fig.lSJ. It 
has been shown in 23] that a simple relation exists between the exponent y and the entropy Cramer function S{x) 
of the Lyapunov exponent of the system (see sec. 0. For the non-universal exponent y the relation y — S'{x) holds, 
where x is the solution of the equation xS'{x) = S{x). Numerical simulations show that in the case of an isotropic, 
short correlated, random velocity, the value of the exponent is approximately y w 3, as shown in Fig.|Sl This implies 



that it is subdominant on the background of the 9 tail coming from the regular dynamics region |0| 



> 



However 



in the general case of a finite-correlated and non-gaussian velocity field, one can imagine a situation where the non- 
universal exponent y becomes smaller than the universal one y < 2, and in this case the 0-angle distribution becomes 
non-universal. 



B. Tumbling time distribution 

For rigid polymers with a fixed size the tumbling time can be defined only through the angular dynamics and 
therefore we will refer only to t^. The evolution time of this system is Tt ~ {s(j>t)^^ so we expect that both the width 
and the maximum of the tumbling time distribution are r^. The tails are related to the probability of passing (clock- 
wise) the angle (j)t after a large amount of independent attempts, and can be estimated as V{t) ~ Q-E{Ds^y^^T 
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FIG. 8; The PDF of the angle 6. 



In |23| it has been shown that the constant E is connected with the ground state energy value of a one-dimensional 
quantum- mechanical system and can be estimated with simple numerical analysis {E ~ 0.45). It is possible to analyze 
the left tail of the tumbling time PDF: 



p{t) cx exp 



2^4(1/2) 



2N-1/3 



(25) 



where K{x) is an elliptical integral of the first kind, and also to give an exact expression for the tail r <^ of the 
PDF. For large values of s/D this tail can be barely observed experimentally or numerically, and the structure of the 
tail strongly depends on the statistics of the chaotic flow, and is therefore non- universal (see Fig.E)). 
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(a) The exponential tail of the PDF. 



(b) The behavior of rt as a function of s'^ D 



FIG. 9: The PDF of the tumbling time t 
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FIG. 10: The behavior of the Lyapunov exponent as a function of D/s. 



V. POLYMER ELONGATION IN CHAOTIC FLOWS 



In this final section we will study elongation statistics of the polymer in the case of random velocity plus mean 
shear. The polymer is not strongly elongated and can be treated as a linear dumbbell, as in sec. lIIII Such a situation 
corresponds to a flexible polymer in a chaotic flow below the coil-stretch transitionlsj- Formally this is the case when 
the maximum Lyapunov exponent is smaller than the inverse relaxation time A < 7, where the Lyapunov exponent 
is the rate of divergence or convergence of two neighboring trajectories. 

The equation governing the system is again eq. where the gradient of the velocity is decomposed into a regular 
shear part and a chaotic part. We can switch to spherical coordinates obtaining an evolution equation for the angles 
(as in eqs. H18(l . ll9|l l and an evolution equation for the modulus of the elongation vector j2M l24l|: 



where in our model we assume {rj{t)rj{t' )) = 2D5{t — t'). Dimensional arguments show that the Lyapunov exponent 
should be proportional to {Ds'^Y^^ (see Fig. IIOII . 

In this situation the polymer spends most of the time in the coiled state, but rare events lead to a large polymer 
stretching, when the flow becomes strong enough. It can be shown, that the right tails of polymer elongation PDF 
have the algebraic form P{R) (x ,23, 32], where the exponent q depends on the parameters of the system. While 

in the thermal noise case the tails were Gaussian, here the probability of observing the polymer stretched is strongly 
enhanced. The algebraic behavior of the tail can be easily explained: the probability of having a local stretching rate 
A > 7 for a large time t decays exponentially with t, and during such events the polymer is stretched by a factor 
which grows exponentially with the time t. 

It has been shown in |23i l3^ that the exponent q is related with the Cramer function of the Lyapunov exponent of 
the flow [23. This function can be found explicitly only in few cases, hence the main aim of this section is to present 
numerical measurements of the Cramer function and of the dependence of the exponent q on the Weissenberg number 
Wix — A/7 below the coil-stretch transition (i.e. WIa < 1). As shown in [2^ the large deviation theory predicts that 
for large averaging times T ^ {Ds'^)^^^'^ = tt the PDF of the Lyapunov exponent is: 



where the function S{x) is the Cramer function. Note that for large but flnite values oi T/tt the measured Cramer 
function Sxix) — — (tt/T) log Pt{x/tt) depends on the time T, but the difference between St{x) and S{x) is 
signiflcant only in the region ^ T/tt 3> 1. Therefore in order to calculate the core of the Cramer function we can 
use the finite time approximation. In our simulations wc considered the ratio T/tt from 10 up to 280. 
In order to connect the Cramer function with the exponent q in the elongation PDF we need to use the Legendre 
transform (see refs. i23|,|3^ for details): 




5t In i? = -7 + - cos^ 6* sin 20 61? 77 , 



(26) 




(27) 



S{x) -xS'{x) +-/ttS'{x) = 
q^S'ix). 



(28) 
(29) 
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FIG. 11: The exponent of the elongation PDF as a function of Wi. 
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FIG. 12: The Cramer function for two different times. 



In and in the behavior of q in the case of zero mean shear has been computed, and the calculations leads to: 

9= I; (7- A), (30) 

where A is the variance of the Lyapunov exponent distribution. In the case of Kraichnan field without mean shear 
the ratio A is proportional to A (see ref. J£| for details), so that q becomes proportional to Wi^^ — 1. In ^34j the case 
of shear turbulence has been analyzed and the scaling is the same as in ea. (|3U|) . 

The exponent q extracted from numerical simulations is plotted in Fig. II II The convergence of the Lyapunov exponent, 
given by the Cramer function, is shown in Fig. 1121 



VI. CONCLUSIONS 



The tumbhng phenomenon [1 171120. l2l| has been studied in the framework of the linear dumbbell model, and some 
universal features of this motion are derived and numerically verified. 

Three different examples of polymers in a linear, steady, plane, shear flow have been studied: i) a flexible polymer 
exp eriencing thermal noise, ii) a rigid polymer in a smooth random velocity gradient above the coil-stretch transition ^ 
|33,|33|), and iii) a flexible polymer in a smooth random velocity gradient below the coil-stretch transition. 
In all three cases, the polymer tumbles aperiodically and the probability density function of the time between two 
successive tumbling events is exponential. While in the case of pure thermal noise, the typical tumbling time is 
determined only by the relaxation process, in the cases of random velocity gradient, the typical rate of divergence of 
two initially close lagrangian trajectories is the most important time scale. 

The analysis of the statistics of orientation leads one to conclude that the phenomenology of these three situations 
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is very similar: in the plane formed by the polymer end-to-end vector and the Z axis there are no relevant scales 
and the polymer spends most of the time in the shear plane. In the shear plane, the dynamics is determined by the 
balance between the shear and the thermal fluctuations. The majority of time is spent nearly aligned to the velocity 
field. Aperiodically the shear induces a tumbling event. 

The elongation dynamics is determined only by the balance between the relaxation and the stretching due to the 
shear flow in the thermal noise case. In the other two cases stretching is determined by the presence of a positive 
Lyapunov exponent. In the case below the coil-stretch transition the tail of the PDF can be determined by measuring 
the statistics of the Lyapunov exponent. 

We would like to thank A. Celani, M. Chertkov, I. Kolokolov and V. Lebedev for many fruitful and inspiring 
discussions and S. Geraschenko and V. Steinberg for many useful comments and interest in our work. We are grateful 
to Chris Carr and Colm Connaughton who helped us to develop the final form of this paper. We acknowledge 
hospitality of CNLS at Los Alamos National Laboratory where this work was partially done. KT was supported by 
RFBR grant 04-02-16520a and INTAS grant Nr 04-83-2922. 



[1] T. T. Perkins, D. E. Smith, and S. Chu, Single Polymer Dynamics in an Elongational Flow, Science 276, 2016, (1997). 
[2] D. E. Smith and S. Chu, Response of Flexible Polymers to a Sudden Elongational Flow, Science 281, 1335 (1998). 
[3] D. E. Smith, H. P. Babcock, and S. Chu, Single polymer dynamics in steady shear flow. Science 283, 1724, (1999). 
[4] J. S. Hur, E. S. G. Shaqfeh, and R. G. Larson, Brownian dynamics simulations of single DNA molecules in shear flow, J. 
RheoL, 44, 713, (2000). 

[5] J. S. Hur, E. Shaqfeh, H. P. Babcock, D. E. Smith, and S. Chu, Dynamics of dilute and semidilute DNA solutions in the 

start-up of shear flow, J. Rheol, 45, 421, (2001). 
[6] S. Gerashchenko, C. Chevallard, and V. Steinberg, Single polymer dynamics: coil-stretch transition in a random flow. 

Submitted to Nature (see also nIin.CD/0404045 I. 
[7] S. Gerashchenko and V. Steinberg , Statistics of tumbling of a single polymer molecule in shear flow. Submitted to Phys. 

Rev. Lett, (see also nlin.CD/0503028 ). 
[8] R. E. Teixeira, H. P. Babcock, E. Shaqfeh and S. Chu, Shear thinning and tumbling dynamics of single polymers in the 

flow-gradient plane, Macromolecules, 38, 581-592, (2005) 
[9] C. M. Schroeder, R. E. Teixeira, E. S. G. Shaqfeh and S. Chu, Dynamics of DNA in the flow-gradient plane of steady shear 
flow: Observations and simulations, Macromolecules, 38, 1967-1978, (2005). 
[10] B. Ladoux, J. Quivy, P. S. Doyle, G.Amouzni and J. Viovy, Direct imaging of singlemolecules: from dynamics of a single 

DNA chain to the study of complex DNA-protein interactions, Science Progress , 84(4), 267, (2001). 
[11] S. Chu, Biology and polymer physics at the single-molecule level , Phil. Trans. R. Soc. Lond. A, 361, 689698 6, (2003). 
[12] S. Manneville, P. Cluzel, J. L. Viovy, D. Chatenay and F. Caron, Evidence for the universal scaling behaviour of a freely 

relaxing DNA molecule, Europhys. Lett., 36, 413, (1996). 
[13] R. G. Larson, T. T. Perkins, D. E. Smith and Chu, S., Hydrodynamics of a DNA molecule in a flow field, Phys. Rev. E, 
55, 1794, (1997). 

[14] Y. Cui and C. Bustamante, Pulling a single chromatin fiber reveals the forces that maintain its higher-order structure, 

Proc. Natl. Acad. Sci., 97, 127, (2000). 
[15] M. Hegner, S. B. Smith, and C. Bustamante, Polymerisation and mechanical properties of single RecA-DNA filaments, 

Proc. Natl. Acad. Sci., 96, 10109 (1999). 
[16] H. Yin, M. D. Wang, K. Svoboda, R. Landick, S. M. Block and J. Gelles, Transcription against an applied force. Science, 

270, 1653 (1995). 

[17] G. J. Wuite, S. B. Smith, M. Young, D. Keller and C. Bustamante, Single molecule studies of the effect of template tension 

on T7 DNA polymerase activity. Nature, 404, 103 (2000). 
[18] A. Groisman and V. Steinberg, Elastic turbulence in curvilinear flows of polymer solutions. New J. Phys. 6, 29 (2004). 
[19] E. J. Hinch and L. G. Leal, The effect of Brownian motion on the rheological properties of a suspension of non-spherical 

particles, J. Fluid Mech., 52, 634, (1972). 
[20] M. Chertkov, I. Kolokolov, V. Lebedev and K. Turitsyn, Tumbling of polymers in random flow with mean shear. Submitted 

to J. Fluid Mech. 

[21] T. W. Liu, Flexible polymer chain dynamics and rheological properties in steady flows, J. Chem. Phys., 90, 5826, (1989). 
[22] R. B. Bird , C. F. Curtiss , R. C. Armstrong and O. Hassager, Dynamics of Polymeric Liquids, Wiley, New York, 1987. 
[23] K. Turitsyn, Polymer dynamics in chaotic flows with strong shear component. Submitted to Phys. Rev. E (see also 
nlin.CD/0501025|. 

[24] Statistics of polymer extension in a random flow with mean shear, M. Chertkov, I. Kolokolov, V. Lebedev and K. Turitsyn, 

submitted to J. Fluid Mech. 
[25] H. C. Ottinger, Stochastic Processes in Polymeric Liquids, Springer, Berlin, (1995). 

[26] A. Celani, A. Puliafito and K. Turitsyn, Polymers in linear shear flow: a numerical study, EuroPhys. Lett., 70, 464, (2005). 
[27] G. Benettin, L. Galgani, A. Giorgilli and J.-M. Strelcyn, Lyapunov exponents for smooth dynamical systems and for 
Hamiltonian systems; a method for computing all of them, Meccanica, 15, 9-30, (1980). 



13 



[28] S.N. Majumdar and C. Sire, Survival Probability of a Gaussian Non-Markovian Process: Application to the T = Dynamics 

of the Ismg Model, Phys. Rev. Lett., 77, 1420 (1996). 
[29] G. Falkovich, C. Gawedzky and M. Vergassola, Particles and fields in fluid turbulence, Rev. Mod. Phys. 73 (4), 913-975 

(2001). 

[30] R. H. Kraichnan, Small-scale structure of a scalar field convected by turbulence, Phys. Fluids 11, 945-963, (1968). 
[31] P. G. De Gennes, Coil- stretch-transition of dilute felxible polymersunder ultra-high velocity gradients, J. Chem Phys, 60 
5030, (1974). 

[32] E. Balkovsky, A. Fouxon and V. Lebedev, Turbulent dynamics of polymer solutions, Phys. Rev. Lett. 84, 4765 (2000); 

Phys. Rev. E 64, 056301 (2001). 
[33] A.Celani, S. Musacchio and D. Vincenzi, Polymer trunsport in random flow, J. Stat. Phys., 118, 529-552, (2005). 
[34] B.Eckhaxdt, J. Kronjager and J. Schumacher, Stretching of polymers in a turbulent environment, Comp. Phys. Comm., 

147, 538, (2002). 



